function y = func_J_m1(q, m, CONSTS)

    k0 = CONSTS.k0;
    a = CONSTS.a;
    p = func_p_from_q(q, CONSTS);
    q1 = func_q1_from_p(p, CONSTS);
    alfa_1 = func_alfa_1(p, CONSTS);
    
    y = (besselj(m+1, k0*a*q1)./((k0*a*q1).*besselj(m, k0*a*q1))) + ...
           ((alfa_1.*m)./((k0*a*q1).^2));
end
